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Abstract 


A common problem for experimental space physicists is the determination 
of the "attitude matrix" T which transforms vectors between representations 
in X and X T coordinate systems according to - T^y* Solution of this 

problem using the standard Eulerian angle techniques is in general tedious 
and circuitous. A straightforward, simple, and efficient solution for the 
transformation matrix is a "double-cross" transformation, such as 

T XX ' “ LA, A x B, Ax (Ax B)j x LA, A x B, Ax (Ax B)J^ f , 

A A 

which is calculated from any two directions A and B which are known in both 
X and X* coordinates. The B direction need be known only well enough to 

A A 

define the plane in which A and B lie. The problem of the intersection of 
two cones as applicable to attitude solutions is also discussed. 



I. 


Introduction 


Experimental space physicists frequently encounter the problem 
of determining the attitude of their sensors on rockets or satellites. 

The purpose of this report is to present a simple, efficient approach 
to the mathematical aspects of the problem of vehicle and sensor 
orientation. 

— ♦ 

It is well known that a (column) vector V , defined in an orthogonal 

X 

coordinate system X f may be transformed to another orthogonal coordinate 
system X having the same origin by a transformation of the form 


V x • T xx' V < l > 

where T is an orthonormal rotation, or attitude, matrix. Due to the 
orthonormal property of T, the inverse transformation is performed by 
the transpose of T: 


V , = T , V 
X r XX X 


, A A 


The problem we shall address is to determine T(A, B) when any two 

A A 

directions, A and B, are known in both X and X representations. (The 

mathematics are, of course, equally applicable in terms of the alternative 

a a 

physical interpretation, that the transformation rotates A and B to new 
positions in a fixed coordinate system.) 

The usual method of expressing T involves Eulerian angles. (See, for 
example, Mechanics texts such as Goldstein [1950], Marion [1965], and others, 
or applications-oriented volumes such as Thomson [1961], Singer [1964] or 
Greensite [1970].) The approach using Eulerian angles is suitable for 
physical problems involving rotations through known angles; however, it 
is difficult to apply in the situation being considered, and is 
cumbersome and expensive for computer applications due to the numerous 
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trigonometric functions involved (29 factors which are either the 
sine or cosine of one of three angles). We shall derive an expression 

A A 

for T(A, B) which involves no trigonometric functions. 

II. The Double-Cross Transformation 

If a matrix M composed of three (or more) column vectors replaces 
V in Equation (1), then T simultaneously transforms each of the vectors 
from one coordinate representation to the other: 

M x = T xx' M x' < 2 > 

Let us use the directions of two vectors A and B, which are known in 
both X and X r coordinates, to construct a set of three orthonormal vectors 

A a 

in each coordinate system. One of several possible sets is A, A x B, and 
a a 

Ax (A x B) . The circumflex denotes a vector normalized to unit length 
— » 

A A A 

A = pH (equivalent to the direction cosines of the vector), and A x B 
denotes that the normalization must be done after the cross product 
operation is performed. The relationship of these vectors is indicated 
in Figure 1. Since the identity of the vectors is independent of the 
coordinate representation, Equation (2) will hold if we take 
M = [A, A x B, A x (A x B)] 

on each side of the equation. Furthermore, M so defined is orthonormal, 
so M ^ = M t . We may immediately solve for T: 

T xx' = ^x ^K' 

r A A A A - r A A A A _t 

= [A, A x B, A x (A x B)] x [A, A x B, A x (A x B)] x , (3) 

Because T can transform any vector between the X and X 1 coordinate 
representations, T determines the relative "attitude" of the two 
coordinate systems. It is clear from Equation (3) that information about 
the magnitude of A or B is unnecessary, and also (see Figure 1) that 
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information about the angle 0^ B has been lost. In applications, therefore, 
because of the sense in which the cross product is taken, A is chosen to 
correspond to the quantity whose direction is more accurately known, and B 
is used only to define the plane in which both vectors lie. Equation (3) 
always has a solution if A and B are not parallel or antiparallel. Obviously, 

the choice of coordinate systems X and X f is completely arbitrary. Even 
in a dynamic system, Equation (3) holds at any instant. The time dependence 
of T may be determined by inserting A(t) and B(t) in the equation. 

Information can be extracted from T concerning the single rotation about 
some appropriate axis that will rotate the X r system into conjunction with 
the X system. The direction cosines of the axis of rotation are the 
components of the eigenvector E of the eigenvalue equation 

(T - [1]) E » 0. 

A 

The rotation angle about the E direction is 

§ = arccos {% (Trace T - 1)}. 

The components of T express the directions of the coordinate axes 
x|, x^j x^ of the X f system in the X system. To illustrate this we take, as a 

A A A A 

special case, A, A x B, and Ax (A x B) parallel to the axes of the X 
coordinate system. Then M^, = [l] and Equation (3) reduces to 
T x = , [A, A x B, A x (A x B)] x = [x{, x£, x-^. 

This special case of the "double-cross" transformation was used by Russell 
(1971), who also appreciated the inconvenience of the Eulerian angle 
approach. 

III. Applications 

In a typical application for rockets or low altitude satellites, A 
and $ might be identified with the directions of the solar and the geomagnetic 
field vectors, X* might be spacecraft coordinates, and X might be geocentric 
equatorial inertial coordinates. 
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A A 

Another choice for A or B might be the total angular momentum vector of 
a spinning spacecraft, since it is often known in both inertial and 
spacecraft coordinates. Definitions of, and transformations between 
common geophysical coordinate systems using the matrix transformation 
approach have been reported in the literature by Russell (1971), or 
may be found in a Goddard Space Flight Center document (1970). 

Not uncommonly, it requires considerable computation to determine 
appropriate vectors in the spacecraft coordinate system. To provide 
a useful example, we shall consider the determination of the solar 

A 

direction, denoted by S , in spacecraft coordinates using three solar 
sensors. Two of these sensors have "look directions" C and D, such 
that the angles 0^ and 9^ between the sensor look directions and the 

A 

solar direction are measured. Thus S may be considered to lie at the 

A 

intersection of two cones with vertices at the origin, with axes C 

A 

and D, and with half-angles 0 and 9 , as depicted in Figure 2. There 

L)b 

A A A 

are in general (for C / D) two solutions for S, and further information 
provided by the third sensor, which must have a look direction not in the 

A A 

plane of C and D, is necessary to resolve the ambiguity. (Note: For 

A A 

rockets or low altitude satellites, C (or D) might, instead, be the 
geomagnetic field vector as determined by an on-board magnetometer. 

In this case, the angle (G^g) between the geomagnetic and solar directions 
need not be measured. This angle is obviously independent of the coordinate 
system and can easily be calculated from ephemeris data) . The problem 
of the intersection of two cones can most easily be solved in the coordinate 
system (X rT ) of Figure 2. The direction cosines of the two solutions 

A 

for S are : 
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where 0 q D = arccos (C • D) , and the ambiguity in S]^ must be resolved 

A A 

with additional information. Equations (4) are inaccurate if C and D are 
nearly parallel or antiparallel. The cones do not intersect if S-^ has 
no real root. To transform S from X ! 1 coordinates to spacecraft (X f ) 
coordinates we use the special case of a "double-cross" transformation: 

M x , = [D x C, C x (D x C), C] x , 

V' -w - M x- 


s x r Vx" S x" M x' M X' t ' ^x" M x' s x" * 

A A A 

At this point S y » may be identified with A or B and the attitude matrix 


T^i found as described previously. 
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FIGURE CAPTIONS 


A A A A 

Figure 1: Three orthonormal vectors A, A x B, and Ax (A x B) created 
from two known vectors A and B. Information about the 
angle is lost, and knowledge of the magnitudes of A or B 

is unnecessary. 

Figure 2: Natural coordinate system for finding the intersection of 

A-f- 

two cones having the same vertex. The two vectors S and 

A. 

S are the solutions sought. Information in addition to 

A A 

D, C, 0^ s and must be provided to distinguish between 

a+ 

S and S . 







